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ABSTRACT 

There has been some controversy over the expression for the so-called 
'interaction energy' due to screening of charged particles in a plasma. Even in 
the relatively simple case of weak screening, first discussed in the context of 
astrophysical plasmas by Salpeter (1954), there is disagreement. In particular, 
Shaviv and Shaviv (1996) have claimed recently that by not considering 
explicitly in his calculation the complete screening cloud, Salpeter obtained a 
result for the interaction energy between two nuclei separated by a distance r 
which in the limit r — > is only 2/3 the correct value. It appears that this claim 
has arisen from a fundamental misconception concerning the dynamics of the 
interaction. We rectify this misconception, and show that Salpeter's formula is 
indeed correct. 

Subject headings: equation of state - nuclear reactions, nucleosynthesis, 
abundances - plasmas - Sun: interior 
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1. INTRODUCTION 

Screening of charged particles in a plasma is important in astrophysics because it 
influences both nuclear reaction rates and the equation of state of the plasma. In general, 
its effect is difficult to calculate. But in the case of classical weak screening of heavy ions, 
in which the Poisson-Boltzmann equation for the electrostatic potential can be linearized, 
the calculation should be relatively straightforward. 

The problem was considered first by Salpeter (1954), who obtained an expression 
for the 'interaction energy' E int (r) between two screened nuclei separated by a distance 
r in a dilute plasma. Recently, Shaviv and Shaviv (1996), in a paper motivated by the 
solar neutrino problem, challenged the procedure by which Salpeter had carried out the 
calculation, and disputed the formula he had derived for E- mt {r) (cf DeWitt, Graboske & 
Cooper 1973, Itoh, Totsuji, & Ichimaru, 1977, Mitler 1977, Alastuey & Jancovici 1987, 
Shoppa et al. 1993). The root of the issue concerns the fact that in evaluating E int (r) 
Salpeter considered explicitly only the contribution to the screening cloud that is associated 
with just one of the screened ions. We explain why this procedure actually leads to the 
correct result in the weak- screening limit, showing how Ei nt (r) can be calculated when the 
entire screening cloud is taken explicitly into account. 

The disputation appears to have arisen in part out of an erroneous view of what the 
interaction energy E- mt should mean. In seeking a resolution it is therefore necessary to 
specify what E int is to be used for. The context in which we discuss the issue is that 
which motivated Salpeter, namely the Coulomb barrier penetration problem associated 
with thermonuclear reactions. This is also the context in which the recent disagreement 
has occurred. In this context, E- mt is the potential energy V appearing in the Schrodinger 
equation describing the relative motion of the interacting ions (nuclei); so from now we 
shall call it V . Misunderstanding is less likely to arise when discussing the equation of 
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state, because in thermodynamics one is forced to differentiate various forms of energy. 

It is common to discuss Coulomb barrier penetration in terms of the JWKB 
approximation to the solution of Schrodinger's equation. In that case, the effect of screening 
is simply to multiply the penetration probability by exp(— V B (Q)/kT), where V B (r) is the 
contribution to V(r) from the screening, k is Boltzmann's constant, and T is temperature. 
Thus it has become common simply to discuss the issue in terms of the value of Vo = V^(0). 
We emphasize, however, that making the JWKB approximation is quite separate from 
making the approximations to the screening energy: the relevance of our discussion is to 
the construction of the Schrodinger equation, irrespective of the approximations that are 
subsequently employed in solving it. 

The Schrodinger equation for the relative motion between two particles is a quantum 
mechanical description of the nonrelativistic balance between kinetic energy and potential 
energy. Kinetic energy is lost or gained, to the benefit or at the expense of potential energy, 
V, which is manifest as a force on the particles. One way of calculating that force is by 
using the principle of virtual work. It is therefore evident that the potential energy V 
appearing in Schrodinger's equation is equal to the virtual work in a realizable motion, as 
Mitler (1977) has pointed out, and it is that quantity whose gradient is the relative force 
between the two screened particles. We demonstrate in this paper how the virtual work is 
correctly calculated. 

After establishing the formulae for the charge distribution and the associated 
electrostatic potential in the weak-screening limit, we argue from a mechanistic viewpoint 
why Salpeter's calculation does yield the correct expression for the virtual work. In the 
subsequent section we reconcile V with the total electrostatic internal energy U of the 
complete two-(screened)-particle system. There follows a discussion of the Shavivs' main 
objections to our argument. 
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Screening is an accumulation of electrons (and light ions) in the vicinity of heavy 
ions, resulting in a systematic deviation of the background charge distribution from that 
uniform state which would have been statistically the most likely had the heavy-ion charge 
distribution been strictly uniform. Consequently, there is a reduction of the entropy of the 
plasma, by an amount S[ nt = (U — V)/T, which precisely accounts for the difference between 
the total electrostatic energy calculated with the entire screening charge distribution and 
the electrostatic interaction energy of each ion (and not its associated screening cloud) 
in the screened potential of the other. The self-energy of the screening cloud which is 
associated with the entropy reduction can be considered as heat, and is not available 
to work at constant temperature on the nuclei to modify their relative kinetic energy. 
Therefore it does not contribute to the potential that appears in the Schrodinger equation 
describing the relative motion of the heavy ions. 

2. WEAK SCREENING 

We consider, as did Salpeter (1954), the screening of ions i, considered as point 
particles with charges Z^e, located at positions r = rj in a nonrelativistic dilute neutral 
plasma at temperature T. Here e is the magnitude of the electron charge. The screening 
cloud is composed of electrons and also, in cases when the screened charge is a (usually 
rare) heavy ion, the relatively light ions of the plasma, whose mean charge number is Z p . 
It is assumed that the motion of the screened ions is slow compared with the motion of 
the screening particles. (For a discussion of 'dynamic screening' of ions for which this 
assumption is not valid, see Carraro et al. 1988.) Then at any instant one may take the 
screening to be the continuous static equilibrium charge distribution in the electrostatic 
field of the ions. Associated with that distribution is an electrostatic potential 0(r) which 
satisfies the Poisson-Boltzmann equation, which may be written 
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V 2 = 47rn c e[exp(+^) - exp(-^)] - 4tt£ Z^(r - r t ) (1) 

in which n e is the equilibrium electron number density, r* is the position of the heavy ion i, 
and S is the Dirac delta function. The reader is referred to Salpeter (1954) for a discussion 
of the conditions under which this equation is valid. 

In the case of weak screening, defined by the condition e<p << kT, equation (1) is 
linearized in e<p/kT. The solution may then be written 0(r) = J2i <f>i, where <pi is what the 
screening potential about the ion % would have been had the ion been isolated from all other 
discrete heavy ions: it is the Debye-Hiickel potential 

4>i{ri\ k) = — exp(-«ri), (2) 

r% 

where r$ = |r — Tj| and k is the reciprocal of the Debye length Ad: 

«m = A5- = ( w £ +z ' ) )' / '- (3) 

In cases where only a fraction a of the ions contribute to the screening, Z p is replaced by 
aZ p in equation (3). The associated contribution to the charge distribution is 

Zen 2 

Pi = — 7 exp(-Kn) + Zie5(r - n), (4) 

the total charge distribution being p = J2iPi- Thus, when linearization is valid, the total 
screening is simply the superposition of the linearized screening of each ion considered 
separately. Indeed, it is sometimes convenient to refer to these two contributions to the 
total screening cloud separately. We shall call them cloud 1 and cloud 2. 
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3. MECHANISTIC DISCUSSION OF THE SCREENED ION 

INTERACTION 

Consider the electrostatic interaction of the screened ions. There are four charge 
components to consider: ion 1, its associated screening cloud, ion 2, and its screening cloud. 
We need to evaluate the repulsive force F12 of ion 1 on ion 2. Thus we consider ion 2 in 
the total screened potential (f)'(ri,r 2 ; k) = <fi — Z 2 e/r 2 of both ion 1 and the entire screening 
cloud. The force is evidently F 12 = —Z 2 eV<f)' evaluated at r = r 2 . However, since the 
electrostatic potential due to the component of the screening cloud associated with ion 2 
is symmetrical about ion 2, it does not contribute to the gradient at the location of ion 
2, and the value of the electric field experienced by ion 2 is the same as the gradient of 
the screened potential of ion 1, which is what Salpeter calculated. Likewise, the value of 
the total potential gradient at r = ri is the same as the value of the gradient of <p 2 . The 
force on ion 2 has magnitude |F i2 | = T\ 2 = — Z 2 ed(p'(ri,r 2 ; K)/dr\ = —Z 2 ed4>\{r\\K)jdr\ 
evaluated at r\ — |r 2 — ri| = r , the distance separating the two ions. The potential energy 
V associated with the interaction is the virtual work: 

V = Z 2 e —— dri = Z 2 e(f) 1 {r; k) = e , (5) 

av\ r 



oc 



in which the integration is carried out at constant temperature, and therefore at constant 

k, and e~ Kr stands for exp(— kt). Of course, one could equally well have considered the 

electrostatic force on ion 1 in the screened potential of ion 2, obtaining the same result. 

The screening contribution V s to the potential energy of interaction, namely the difference 

between V and the potential energy Z\Z 2 e 2 /r of the corresponding unscreened interaction, 

is 

Vs(r) = ^V--l), (6) 

r 

whose value at r = is Vq = —ZiZ 2 e 2 K. This formula is equivalent to that obtained by 
Salpeter (1954). 



Of course there is also an attractive force between the screened ion 1 and the screening 
cloud about ion 2, whose magnitude is J- c \. But since in the linearized approximation 
the screening cloud remains distributed symmetrically about its ion, none of that force is 
transmitted to ion 2. It is transmitted instead to the rest of the plasma. Therefore, in this 
method of calculation, the force of both ion 1 and its screening cloud on the screening cloud 
of ion 2 does not contribute to V. We note that the forces exerted by the plasma on the 
two components of the screening cloud must be equal and opposite, because no net force 
can be exerted by the otherwise uniform plasma on the screened ion pair. 

In reality, the different forces on the ion and its associated screening cloud polarize 
the system, causing the ion and its screening cloud to interact electrostatically, thereby 
reducing somewhat the magnitude of V . But that reduction is nonlinear in e<j)/kT, and is 
formally negligible in the weak-screening limit. 

It is interesting to note that if it were to be assumed that the polarization is simply 
a small displacement of the screening cloud relative to its ion, without distortion, such 
that the effect is simply to transmit the force on the cloud to the ion, V would become 
Z\Z2e 2 r~ l (l — |/tr)e _Kr , and Vq would be — ^ZiZ 2 e 2 K, which is the value preferred by 
Shaviv and Shaviv (1996). However, such an assumption is inconsistent: if the ions are 
labeled such that Z\ > Z 2 there is no possible nontrivial displacement of the undistorted 
screening cloud of ion 2 that can transmit to that ion a force equal to T c \. Any consistent 
generalization of the weak-screening formula must take due account of the distortion of 
each screening cloud by the electrostatic attraction of the neighbouring screened ion. 
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4. THERMODYNAMICAL DISCUSSION OF THE SCREENED ION 

INTERACTION 

We now regard the screened ion pair discussed in the previous section as an integral 
part of the plasma. Our starting point is the electrostatic energy of their interaction, which 
contributes to the internal energy of the plasma. It is given by 

[/=l|0(r)p(r)d 3 r, (7) 

the integral being over all space. The components <pi and pi of and p are given by 
equations (2) and (4). Evaluation of the integral is straightforward, and is outlined in the 
Appendix. The result is 

U(r,T) = U- mt (r,T) +U oa = M^(l - l - K r)e~^ + */«,, (8) 

r 2 

in which r is again the separation of the ions, and U 00 {T') is independent of r. The 
contribution to U from the screening is U s — U — Z\Zi^ jr. 

In order to evaluate the potential V, we apply the first law of thermodynamics: 

dU = TdS - F 12 dr, (9) 

where S is the entropy of the screened ion pair; it is given by 

where F = U — TS is the Helmholtz free energy. The second of equations (10) integrates at 

constant r to 

r°° U(r T') 

F(r,T)=Tj T ^^dT'. (11) 

From equation (9) and the definition of F one obtains dF = —SdT — jF 12 dr, from which it 
follows that 
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The force between the two screened particles is minus the isothermal gradient of the free 
energy, irrespective of the degree of screening, as DeWitt, Graboske and Cooper (1973) 
and Mitler (1977) have pointed out. It is not minus the isothermal gradient of the internal 
energy. Thus the potential V in Schrodinger's equation is always F int = F(r,T) — F(oo,T), 
whether the screening is weak, intermediate or strong. With this identification, substituting 
the weak-screening expression (8) for U into equation (11) yields V = F int = ZiZ 2 e 2 r~ l e~ KT , 
which recovers the final expression in equation (5), and hence Salpeter's expression (our 
equation 6) for V s . The repulsive force on ion 2 in the electrostatic potential of both ion 1 
and the (entire) screening cloud, according to equation (12), is given by 

Fi2 = ^f-(1 + Kr)e- Kr , (13) 

and the entropy, which can be calculated from equations (10), is 

S(r,T) = S 00 (T)-^^e-, (14) 

where S^ is the limit of the entropy at temperature T as the separation between the ions 
tends to infinity. 



5. DISCUSSION 

In the weak-screening limit, the Poisson-Boltzmann equation (1) is linearized, so that 
the screening cloud surrounding two stationary (or slowly moving) heavy ions can be 
represented by a superposition of the two spherically symmetrical screening clouds that 
would accumulate independently about each ion considered in isolation. The concentration 
of screening particles in the potential well of the two ions reduces the electrostatic energy 
of the system; it signifies an increase in order in the plasma, and therefore a reduction in 
entropy. 
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When considering the consequences of screening to barrier penetration, it is necessary 
to calculate the potential energy V of the mechanical interaction between the ions, whose 
derivative with respect to the ion separation r is the force between those ions. That force is 
purely electrostatic: the force on ion 2 is simply the product of the charge of ion 2 and the 
electric field due to both ion 1 and the total screening cloud. However, since the electric 
field of the component of the screening cloud associated with ion 2 vanishes at its centre of 
symmetry, where ion 2 is located, the value of the total electric field at this location, and at 
this location only, is the same as the value of the field due only to ion 1 and its associated 
screening cloud. The electrostatic force between the ions can be calculated either as the 
isothermal derivative of the free energy, or as the product of the charge of ion 2 and that 
part of the screened charge of ion 1 that is contained in the sphere of radius r about ion 1 
(obtained by integrating p, given by equation 4, over the volume of the sphere) divided by 
r 2 , or the equivalent expression obtained by interchanging the labels 1 and 2. The potential 
from which that force is derived is the potential V that appears in Schrodinger's equation; 
it is the same as that adopted by Salpeter. 

In his referee's report on the preceding sections of this paper, G. Shaviv, in consultation 
with N. Shaviv, has claimed that our argument is fallacious. We address here just his 
principal two objections, in the hope of illucidating the physics involved. 

The essence of the first objection is based on the fact that, in what Shaviv calls 
the adiabatic picture, the total screening energy, namely the difference between the total 
interaction energy U- mt given by equation (A3) and the interaction energy Z\Z2e 2 /r of 
the bare ions, is (at zero ion separation) 3/2 of the value of K(0) given by equation (6). 
Shaviv complains that he cannot understand why that should not imply that our value of 
Vq = V s (0) must be multiplied by the factor 3/2. 

The reason why it must not can be understood as follows: The difference between U- int 
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and V is the interaction energy U c \ between the two cloud components, associated with 
which is a contribution to the entropy of S c \ = T~ l U c \ of the screened-ion system. So far as 
the mechanics of the ions is concerned, that cloud-cloud interaction energy is essentially 
heat, and it is not available to contribute to the force between the ions in the isothermal 
motion. The plasma exerts a force on each screening cloud which balances the electrostatic 
force exerted on the cloud by the other screened ion. The energy U c \ is exchanged with 
the background plasma, which acts as a heat bath. It is the remaining energy, U- Ult — U c \, 
namely the interaction free energy Fi nt , whose isothermal derivative with respect to r is the 
mutual force between the ions, and which must therefore be identified with the potential V 
in Schrodinger's equation. 

We suspect that part of the confusion experienced by Shaviv has arisen from an 
incorrect assumption that the actual motion of the screened system is adiabatic. To be sure, 
one can calculate the force between the ions by considering a virtual adiabatic displacement. 
The force is evidently the derivative with respect to r of the total electrostatic interaction 
internal energy U- m t at constant entropy S, because in such a virtual displacement there 
is no heat exchange with the plasma. The outcome is, of course, the same as that given 
by equation (12), since (dU/dr)s = (dF/dr)?- This may be verified by differentiating 
expression (A3) for U- int at constant S (which is given by equation 14), yielding equation 
(13). In a virtual displacement corresponding to that derivative, temperature varies, so 
the displacement does not correspond to the actual motion. The real motion is isothermal, 
because the screened heavy ions move slowly in the heat bath provided by the rapidly 
moving electrons and light ions of the surrounding isothermal plasma. Therefore, in 
computing from this result the potential V, the integral / T\idr must be carried out at 
constant temperature, which yields equation (5). 

We observe that the result of Shaviv and Shaviv (1996) is obtained by replacing the 
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force with the isothermal derivative of U- int . That is evidently an incorrect procedure, 
because it ignores the heat exchange with the plasma. Moreover, it would predict that the 
force between the screened ions is attractive when nr > 1 + \^3. It is most straightforward 
to see that that cannot be so by considering the special case Z\ = Z 2 ; then the effective 
charges of each screened ion experienced by the other are evidently identical, and the force 
between them must therefore necessarily always be repulsive. 

Shaviv's second objection rests on the mistaken belief that in our calculation of the 
force exerted by ion 1 on ion 2 we do not take into account the electrostatic potential of the 
screening cloud associated with ion 2, and that consequently that force is not necessarily 
equal and opposite to the force exerted by ion 2 on ion 1. He points out that the force on a 
bare ion 2 in the electrostatic field of the screened ion 1 is different from the force on ion 1 
in the electrostatic field of the bare ion 2, which, to the best of his knowledge, contradicts 
Newton's third law. 

We do not take issue with the second of these statements, because it does not apply 
to our calculation. But we do take issue with the first. That we do not ignore the 
contribution to the electrostatic potential from the screening cloud associated with ion 

2 is most readily seen in the Appendix, for its contribution to the interaction energy 
is explicitly the third term of the integrand in equation (Al). We point out in section 

3 why the electrostatic potential of screening cloud 2 does not contribute to the force 
on ion 2: we did take it into account, but it contributed nothing. And likewise, the 
potential of cloud 1 does not contribute to the force on ion 1. Contrary to Shaviv's 
claim, the situation described by us is symmetric with respect to the two ions, and action 
and reaction are equal and opposite. Indeed, Fi 2 = —Z 2 ed(f)i(ri;K)/dri\ ri=r dr/dr 2 and 
F 2 i = — Zied(f> 2 (r 2 ; K)/dr 2 \r 2 =rdr/dri, where, as previously, r = |r 2 — r x |; and since 
Z\e(j) 2 (r; k) = Z 2 e<f)i(r; k) and dr/dri = —dr/dr 2 it follows that F 2 i = — F 12 , in agreement 
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with Newton's third law. It is straightforward to demonstrate that the force of the plasma 
on the screening cloud of ion 2 is equal and opposite to the force of the plasma on the 
screening cloud of ion 1, confirming that the plasma exerts no net force on the screened ion 
pair. 

Finally, one must bear in mind that although the principles of our discussion are valid 
quite generally for essentially static screening, the detailed formulae are valid only in the 
limit of weak screening. In reality, each ion together with its screening cloud are polarized 
by the presence of the other, which modifies the interaction free energy. Moreover, the 
polarization depends not only on the separation of the ions, but also on their velocities. 
To take due account of the polarization, a proper quantum mechanical description of the 
electron distribution is required. That is beyond the scope of the present discussion. 

6. CONCLUSION 

In the weak-screening limit, Salpeter's expression, essentially equation (5), for the 
interaction energy that appears in the Schrodinger equation governing the relative motion 
of two screened ions is correct. In particular, the contribution from screening at zero 
separation is Vq = —ZiZ 2 e 2 /\n, where Z\e and Z 2 e are the ionic charges and Ad is the 
Debye length of the light component of the plasma. The more recent incorrect suggestion 
that Salpeter's value for Vq should be multiplied by a factor 3/2 seems to have arisen from 
the false assumption that the electrostatic force on a screening cloud due to the other 
screened ion, evaluated in the linearized approximation, is transmitted directly to the ion 
it surrounds, rather than to the surrounding plasma, and from a failure to appreciate that 
heat is exchanged between the screening cloud and the surrounding isothermal plasma as 
the distance r between the ions varies. When the heat exchange is correctly taken into 
account, it transpires that it is Helmholtz free energy, and not internal energy, that must 
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be identified with the potential energy V in Schrodinger's equation. Not only do the free 
energy and the internal energy of the screening differ in magnitude at zero ion separation - 
-^int = f^int at r = - but they also have different functional dependences on r. Therefore, 
to use the interaction internal energy U- m t rather than the interaction free energy i*j nt for 
the potential V in Schrodinger's equation leads to yet further errors when one includes 
corrections from the Liouville-Green expansion about the JWKB approximation to the 
barrier penetration probability. 

We are grateful to D. Lynden-Bell for illuminating discussions. We thank our first 
referee, G. Shaviv, for insisting that our discussion in sections 1-4 does not make the case, 
thereby causing us to add section 5 which we trust illucidates the matter satisfactorily. We 
also thank H. DeWitt, the second referee, for stating unambiguously that our discussion is 
correct. 



A. APPENDIX 

After making the substitutions <f) = <pi + <p2 an d p — P\ + P2 in equation (7), and 
using equations (2) and (4) for fa and pi, one can identify two distinct groups of terms: 
those representing the sum Uoo of the electrostatic energies of each of the two screened 
ions considered in isolation, and those representing the interaction energy between the two. 
The former does not vary with the separation r between the ions, and therefore makes no 
contribution to either the force between the ions or the free-energy excess F- mt over the free 
energy at infinite ion separation. Therefore we need consider only the latter group of terms. 
It is given by 

1 „ „ o f f 2k' 



47rr 1 r 2 



p "" 1 p '"» A 

S(r - r 2 ) H S(r - ri) 

r\ r 2 



d% (Al) 
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where r\ and r 2 are the distances of the volume element of integration from the respective 
ions, and the integration is over all space. To evaluate the first term, U c \, which represents 
the excess of the electrostatic self-energy of the complete screening cloud over the value 
when the ions are infinitely separated, it is convenient to adopt ellipsoidal coordinates 
(£,?7,0), where now is the azimuthal angle about the axis upon which the ions lie, and 
r £ — r i + r 2) r V — r i — r 2- The appropriate part of the right-hand side of equation (Al) 
then reduces to 

I fl f-OO I 

U d = --Z l Z 2 e 2 n 2 r / dri / e" Kr5 d£ = — Z 1 Z 2 e 2 Ke~ Kr . (A2) 

4 J-i J\ 2 

The remaining contribution is the sum of the electrostatic energies of each of the two ions 
in the screened potential of the other. Its value is evidently Z 1 Z 2 e 2 r _1 e _Kr . Hence 

U int = Z 1 Z 2 e 2 r~ 1 (l - l -Kr)e~ Kr . (A3) 

For completeness, we record that U^ = —\{Z\ + Z 2 )e 2 K. The total electrostatic internal 
energy U is U- m \, + f/oo, as stated by equation (8). 
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